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An experimental method is proposed which allows the complete determina- 
tion of the complex reflection coefficient for absorptive media for positive and 
negative values of the momenta. It makes use of magnetic reference layers 
and is a modification of a recently proposed technique for phase determina- 
tion based on polarization measurements. The complex reflection coefficient 
resulting from a simulated application of the method is used for a reconstruc- 
tion of the scattering density profiles of absorptive non-magnetic media by 
inversion. 



I. INTRODUCTION 

Neutron specular reflection has become a common tool in materials sciences 0. Never- 
theless, there are still severe difficulties in the interpretation of reflectometer experiments 
because of the so-called phase problem. This problem refers to the incompleteness of the 
data obtained in standard experiments in structure physics where only the intensities 
of the reflected waves are measured but not the corresponding phases. In specular reflection 
both the modulus and the phase of the reflection coefficient are needed for an unambiguous 
reconstruction of the surface profiles @||. 
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Several solutions of this so-called phase problem have been proposed for neutron specular 
reflection PJ6|-|T0|| . The recently proposed reference layer methods based on polarization 



measurements |p|,[10[] are of particular interest because they also work in the total reflection 
regime and allow the unique reconstruction of surface profiles of magnetic samples. Despite 
the great effort in theory and modelling, only the phase determination by the method of 
Majkrzak and Berk [0 was tested experimentally JXX[] for a specific case. 
The available methods of phase determination [e.g. [f7t^,|ITi|l2|1 for non-absorptive as well as 
for absorptive samples provide us with the full (modulus and phase) reflection coefficient R{q) 
for positive values of the wave number q perpendicular to the surface of the sample. On the 
other hand, the reconstruction of the surface profiles by inversion requires the knowledge of 
R{q) over the positive and negative range of q- values. This is no problem for non-absorptive 
samples because in this case the reflection coefficient R{q) satisfies 

R(-q) = rf(q) , (1) 

and R(q) at negative g-values is directly given in terms of R(q) at positive q. However, in 
the presence of absorption relation (JIJ) does not hold, and therefore a reconstruction of the 
surface profiles is not possible without further input. To our knowledge, so far no procedure 
has been proposed for providing the missing information. 

In this work we propose a modification of the reference layer method for non-magnetic sam- 
ples based on polarization measurements ||. The modification consists in the interchange 
of the positions of the sample and the reference layer, which allows the determination of the 
reflection coefficient for positive and negative (/-values. The method works also in the total 
reflection regime (in contrast to [f7|,^|,[lT|); it is only limited by the sample thickness because 
of the corresponding Kiessig oscillations which must be resolved. 

In section II we present the method and derive the basic relations for the phase determina- 
tion. A realistic example is given in section III, where we demonstrate by simulations the 
feasibility of the method. The problem of thick samples is also discussed. A brief summary 
and concluding remarks are given in section IV. 



II. THE METHOD 



We consider the arrangement of an unknown non-magnetic sample and a magnetic reference 
layer mounted on a substrate (e.g. a Si-wafer) as shown in Fig. 1. We assume that within 
the reference layer (0 < x < a) there is a magnetic field B(x) aligned with the surface of 
the sample. The direction of this field is taken as the z-axis, which is chosen as the axis 
of spin quantization. The direction of propagation perpendicular to the surface defines the 
x-axis in a right-handed coordinate system. Such a magnetic reference layer may consist of 
a ferromagnetic stratum, e.g. a Fe-, Co- or Ni-layer (cf. Fig. 1). We assume that the sample 
as well as the substrate are field-free. However, the further considerations remain still valid 
as long as the condition |B(x)| = for x — > oo is satisfied. 

The arrangement of Fig. 1 differs from that of §J only by the interchange of sample and 
reference layer and therefore we can make use of the same relationships for the description 
of the reflection. The position of the sample on top of the reference layer and the substrate, 
however, leads to a completely new situation as compared to previous proposals for the 
solution of the phase problem. This is easily seen from the expression for the reflection 
coefficient R tot of the whole arrangement 



E(q)p L (q) + R L (q) 
1 - R R (q)p L {q) 

Here the reflection properties of the sample enter in terms of the reflection coefficients Rl 
and R R and the quantity 

E(q) = T L (q)T R (q) - R L (q)R R (q) , (3) 

where Tl, Tr are the transmission coefficients of the sample. The quantity pi is the reflection 
coefficient of the reference layer plus the substrate and the indices L, R refer to incidence 
from left and right, respectively. Because of the rather intricate dependence of R 1 ^ on 
the reflection properties of the sample the separation of its unknown reflection coefficient 
becomes more involved. 
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Following the procedure of || we introduce the quantity 



Dtot 

s= !f' (4) 



which can be expressed in the form 



' P^ + ^l + po' W 

where P = (P x , P y , P z ) and P° = (P°, P°, P°) are the polarizations of the reflected and the 
incident beam, respectively. The indices ± refer to neutron beams polarized parallel (+) 
and antiparallel (— ) to the magnetic field. All quantities on the right-hand side of Eq. fl5|) 
are measurable and therefore s can be determined experimentally. 
Use of Eq. (pD for P±* in Eq. (||) leads after simple algebraics to the relation 

_ Rh{q) = (pt ~ sp~ L ) - R R pjpl(l - s) , , 

' E(q) {l-s)- R R (p L - apt) ' 1 ) 

The term on the left-hand side of Eq. (||) is equal to Rx(—q), i.e. to the right reflection 
coefficient at negative q- values (cf. |IByrjj). This is easily seen from the Jost solutions of 
the bare sample. Making use of the linear independence of the Jost solutions fi,,R(q, x) with 
asymptotic forms 

, x \ \e iqx + R L {q)e- tqx ] /T L {q) fori^-oo 
h(q,x) = l (7) 

e iqx for x — > +oo 

and 

{e~ iqx for x — > — oo 

(8) 
[e-' iqx + R R (q)e iqx ] /T R {q) for x -> +oo 

one finds 

/«(-?, = T L (q)f L (q, x) - R L (q)fR(q, x) . (9) 

Comparing Eq. (^]) in the limit x — > oo with the corresponding asymptotic form (^) with g 
replaced by — g we obtain 



and therefore, from Eq. (P). 



(1-s) -R R {q){p L -spl) 

This is a relation involving measurable quantities which expresses the right reflection coef- 
ficient at — q (q > 0) in terms of its value at q. 

For a full determination of the reflection coefficient Rn{q) at positive and negative q- values 
(i.e. for obtaining two independent equations for the two unknowns Rn{q) and Rni—q)) at 
least two sets of measurements with different reference layers are necessary. Denoting the 
two measurements with the upper indices w and ^ one obtains from Eq. fllTD the quadratic 
equation for Rn(q) 

a + (3R R (q)+ 1 R 2 R (q) = (12) 

with 

a = (1 - s^)(p^ - 8 MpZ {a) ) " (1 " * (a) )(p£ (6) - s (6) PZ (6) ) , (13) 

P = (Pl ib) - s^ P + L {b) )(p + L ia) - » (a W o) ) " (Pl ib) ~ s (b) Pl (b) ){pl {a) - * ( "VJ (a) ) 
+ (1 - S W)(1 - s^)pt {a) pl {a) - (1 - - s^)pt {b) p- L {b) , (14) 

7 = (1 - S (a) )(p L W " s^p + L ^)p + L (a) p- L {a) - (1 - , (6) )(pZ (a) - ^ (a) )p^ ( V fe) • (15) 
The quadratic equation ( |12|) has two roots 

Similarly to the procedure of Ref. || the physical solution Rn(q) can be selected either 
by continuity requiring that the phase $>(q) satisfies <j)(q = 0) = — 7r or by the condition 
tr = \Rr\ 2 < 1. The latter may fail at some q- values and must be accompanied by continuity 



in certain momentum regions. The selection by continuity implies measurements of s over 
the whole range of momenta. The solution Rji(—q) at negative q- values is given by Eq. (|TTD. 
We note that the method relies solely on measurements of the polarization of the reflected 
beam. In particular, even for absorptive samples, no transmission measurements are needed 
for the determination of R{—q) as might be conjectured from the general formalism 



III. EXAMPLE 

We test the method by a simulation using the realistic example shown in Fig. 1. The sample 
consists of a 30 nm Au- and a 10 nm Cd-layer, where the latter is strongly absorptive. To 
determine the reflection coefficient for positive and negative g-values we consider measure- 
ments with two different reference layers. Both reference layers are composed of a 15 nm 
Cr-layer and a 15 nm thick film of ferromagnetic material magnetizied up to saturation and 
mounted on a Si-wafer. The required difference in the reflection properties of the reference 
layers can be achieved by using Fe in measurement (a) and Co in measurement (b) for the 
ferromagnetic layer. The magnetization of this layer will generate a magnetic induction also 
outside the ferromagnetic film which we assume, however, to be small enough not to affect 
the neutron beam. For simplicity we set it equal to zero. 

The reflectivity r^°| = |-R^+| 2 and the polarization components P x ,P y ,P z for the arrange- 
ment of Fig. 1 have been calculated for both reference layers at all positive g-values. To 
illustrate the typical observables which have to be measured for the phase determination 
we show in Fig. 2 the reflectivity r l £\_ = |-R^+| 2 and the components of the polarization 
P x ,P y , P z for the reference layer (a) assuming a fully polarized incident beam in x-direction. 
It is evident that in the given example measurements in the range < q < 1.0 nm" 1 are 
certainly feasible. 

Following the procedure outlined in section II we obtain two solutions for Rn(q) which are 
displayed in Fig. 3. As expected, in a wide range of positive g-values only one solution 
satisfies \Rn(q)\ < 1 and is therefore physically admissible. Applying continuity the region 
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of uniqueness can be extended over all positive g-values. In addition one obtains via Eq. (|TTD 
unique values of Rft(-q). Hence the measurements of the polarization at positive q- values 
is sufficient for the determination of the complete reflection coefficient on the whole g-axis. 
To extract the surface profiles from the extracted reflection coefficients Rn{q) we apply 
the Marchenko inversion procedure used in Ref. [IB|, where it was shown to yield correct 
reconstructions of the complex input potential, although a rigorous mathematical proof of 
this is lacking. The Marchenko inversion procedure is usually formulated for Ri(q), while 
our procedure provides Rii(q). Because of the specific choice of our coordinate system we 
can also use Rn{q) in the formulation of Ref. |l3j but we have to change the argument 
of the evaluated potential from x to —x. This procedure has been applied to analyse the 
numerically extracted values of R^q). The results for the scattering length density profile 
are displayed in Fig. 4. It is obvious that the method reproduces the original potential 
within the resolution determined by the maximum available g-value. 

Since the method is strongly related to the reference layer method presented previously || 
its stability with respect to experimental uncertainties, e.g. roughness of interfaces and 
measurement errors, is similar to that of method || . This is also true for the reconstruction 
of the surface profiles. For more details on the stability we refer to our previous studies (cf. 
ill)- 



IV. CONCLUSIONS 

We have proposed a method involving a magnetic reference layer which, for the first time, 
allows the determination of the full complex reflection coefficient at positive and negative 
g-values. Thus a reconstruction of surface profiles of absorptive non-magnetic samples be- 
comes feasible. The novel procedure is a modification of the previously proposed method 
||. The main difference is the interchange of the position of the sample and the reference 
layer. Because of this change the reflection properties of the sample enter in a more in- 
volved way and require two sets of measurements using two different reference layers. The 
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additional effort is remunerated by providing us with an additional relationship from which 
the reflection coefficient at negative g-values can be obtained. This novel feature allows the 
determination of absorptive surface profiles by inversion, - an open problem up to now. In 
a numerical example it is demonstrated that the reconstruction of real and imaginary parts 
of the scattering density profiles is feasible. 

As in the previously proposed reference layer methods [0-|9| the thickness of the layers 
is reflected in the g-dependence of the measured data. Specifically, the so-called Kiessig 
oscillations must be resolved in experiment. Hence, the degree of the monochromaticity 
as well as the resolution of the diffraction angle of the incident beam set a limit on the 
thickness of the studied system. The roughness of the interfaces has similar effects and 
results in additional uncertainties in the deduced reflection coefficient (cf. 0). 
In summary the phase problem of neutron reflection has been solved for absorptive non- 
magnetic samples. The method seems experimentally feasible for sample thicknesses up to 
100 nm and complements the reference layer methods |9||ll| recently presented. 
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FIGURE CAPTIONS 



Figure 1 

Experimental arrangement for measuring the complex reflection coefficient Top: Arrange- 
ment of the layers. Bottom: The potential profile; the real part is represented by the solid 
lines, the imaginary part by the dashed area. The dotted lines represent the effective poten- 
tials experienced by neutron beams polarized parallel and antiparallel to the magnetic field 
B. 

Figure 2 

Simulated reflectivity and polarization data for the arrangement of Fig. 1 with a magnetized 
Fe-film in the reference layer. The incident beam is assumed to be fully polarized in the 
x-direction. The reflectivity r^j. = |-Ri°|| 2 (a) and the polarization components P x (b), P y 
(c) and P z (d) of the reflected beam are shown. 

Figure 3 

(1 2) 

The roots R R ' of the quadratic equation as obtained from the measurements of the polar- 
ization components. Top: the reflectivities r^' 2 ** = li?^' 2 "*! 2 . Bottom: the absolut reflection 
phases <j) R ' are shown assuming 0^ = — tx at q = 0. The solid curves correspond to the 
physical solution. 

Figure 4 

The potential profile of the sample obtained by inversion of the extracted values of the 
reflection coefficient in the momentum range —1.0 < q < 1.0 mil" 1 shown in Fig. 3. The 
real and the imaginary parts of the reconstructed profile are shown by thick solid and dashed 
lines, respectively. For reference the original profile is shown by thin solid lines. 
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